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1 Introduction 

The folklore reading "Every symmetric matrix can be brought into diagonal 
form by a suitable rotation" is strictly valid in the positive definite case only. 
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In the Lorentz signature case, however, one needs additional assumptions to 
get this result. It is generally believed that these assumptions do not rep- 
resent a real restriction, and this is justified, e.g., for the energy-momentum 
tensor: all physically sensible form of matter can be written with an energy- 
momentum tensor in diagonal form. 

It is the purpose of the present paper to show, that, nevertheless, impor- 
tant examples exist where this folklore-statement leads to incorrect results. 
Of course, some examples of this kind already exist. Most notably, the Kerr 
metric, and all its generalizations, cannot be brought into diagonal form in 
holonomic coordinates due to the fact that its timelike Killing field fails to 
be hypersurface-orthogonal. However, there is a widespread feeling that, in 
the case of the metric depending only on one coordinate, this diagonalization 
can always be achieved. Here we give a class of solutions where this rough 
diagonalization is not possible. This is essentially important if one wishes to 
find all solutions of a prescribed symmetry type. In this respect, this is a 
continuation of and [Tj. 

The paper is organized as follows: Section 2 presents the deduction of the 
Bianchi type I vacuum solutions of Einstein's gravitational field equations in 
such a detailed form that it becomes clear, why the Kasner solution [3], 
here cited from [5J,^ is really possible in diagonal form.^ Section 3 gives the 
analogous calculation as section 2, but now with a changed signature. In 
the final section 4 we discuss the case of a metric with only one off-diagonal 
term. We will show that the vacuum Einstein equations for this kind of 
metric, that is a system of four nonlinear differential equations with four 
unknown functions, can always be reduced to a system two equations with 

^According to the Kasner solution published in 1921 was already known to Levi- 

Civita in 1917, cf. also 

■^To the question when symmetric matrices cannot be diagonalized see 
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two unknown functions. Moreover, in some cases we will be able to further 
reduce the system to only one equation in one unknown function. We will 
also show one explicit solution of the system, that is a vacuum space-time 
which cannot be diagonalized. 



2 The Kasner solution 



The general metric for a Bianchi type I model reads 



(1) 



where Qapix) dx'^dx^ is the positive definite spatial metric. We want to find 
out all vacuum solutions of the Einstein field equation of the form of met- 
ric (Q). The result reads, cf. section 11.3. of [S]: 



ds^ = -de + t'^dx' + t"^dy^ + e-dz" 



(2) 



with 



p + g-hr =p2 + g2 + r2 e {0, 1} . (3) 

Proof:^ This is a well-posed Cauchy problem. We take \t = 0] as initial 
space-like hypersurface. It is a 3-flat space. Therefore, we can take without 
loss of generality 

/ 1 \ 



(4) 



1 

Vo I) 

otherwise a coordinate transformation involving only the 3 spatial coordi- 



z suffices to reach that. 



nates x^ = x, x"^ = y, x^ 

''it is a fully standard proof, but here we need the details to see the distinction to the 
other cases to be discussed below. 
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The second fundamental tensor at [t = 0] is ka/3{0), where 

kaP = ^ gal3 (5) 

represents a symmetric matrix. The remaining freedom of spatial coordinate 
transformations keeping valid the equation is just the orthogonal group 
0(3), see Appendix A for the proof that the 3 parameters of 0(3) suffice^ to 
reach 

A;i2(0) = fci3(0) = A:23(0) = 0. (6) 

Therefore, due to the compactness of the rotation group we can assume 
without loss of generality that the second fundamental tensor at [t = 0] has 
diagonal form. 

The vacuum Einstein field equations are equivalent to Rij = 0, where Rij 
represents the Ricci tensor, and the index of the coordinates covers the 
values i = 0,1, 2, 3, and t = Then it turns out that the 3 equations 

Ri2 = Ri3 = R23 = (7) 

suffice to maintain the diagonal form of the metric Qapit) for all times. 

Up to now we have shown that all vacuum metrics of the form can be 
written as 

ds^ = -de + e^^Wrfa;^ + e^P^'^dy^ + e^^^rf^^ _ ^g) 

Inserting this metric into Rij = it turns out that up to trivial rescalings, 
the one-parameter set of solutions defined by eqs. (j2l El) cover the set of all 
solutions. 

"^In fact, we have proven even more: already the connected component of the unity 
element of 0(3), namely the SO{i), is enough to get that result, but we do not need this 
additional property here. 
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In the final step one observes, that —dt'^ + dx"^ and —dt'^ + t'^dx'^ are both 
locally flat, and therefore, one may omit the case p = q = r = from eq. 
without loosing any solutions. 

Result: Every cosmological Bianchi type I solution of the Einstein vacuum 
field equations, i.e., every solution of the form can be written as (j21) with^ 



3 The signature changed Kasner solution 

In this section we want to deduce the consequences of another signature in 
the metric. First of all, one would be tempted to go just the same way as 
before: Looking at eqs. (j21 El), one can perform an imaginary rotation of x 
and y. After rewriting ds"^ as —ds"^ and renaming the coordinates one gets: 



It holds: If eq. is fulfilled, then metric (fTI]|) represents a vacuum solution 
of the Einstein field equations. Contrary to the positive definite case, a 
permutation between p, q, and r is no more generally possible. It remains 
only the permutation between q and r. So, the set of solutions in metric ()lUp 
can be parametrized by eqs. (^01 HH |121) with < < vr. 

We will now carefully look for the question whether these solutions rep- 
resent all solutions of the form of metric (QJ). 

3.1 The diagonal ansatz 

The diagonal ansatz analogous to eq. (jH)) reads 



p + g + r=p^ + g^ + r^ = l . 



(9) 



ds^ = dx^ + x^^'dy^ + x^'^dz^ - x^Pdt^ 



(10) 



ds' = dx' + e^-^^^^dy' + e^^^^^dz" - e^^^^^^dt^ . 



(11) 



geometric parametrization of this set is given in Appendix B. 
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Inserting this metric into the equation Rij = 0, one gets as expected, again 
just the known solutions (fTUj) with Q. 

3.2 The non-diagonal ansatz 

We are now only interested to show that truly non-diagonal metrics really 
exist, and we do not intend to exhaust all of them in the present paper. 
Therefore, we restrict to those metrics, where only one off-diagonal element 
of Qapix) is different from zero. 

An off-diagonal component between two space-like directions can be made 
vanish by the procedure shown in Appendix A. So, this essential off-diagonal 
component must exist between one space-like and one time-like direction. 
This leads us to the following ansatz for the metric: 



where A{x) > and gmn{x) is negative definite. We count the coordinates 
= z and = t, so the indices m, n run from 3 to 4. We write the 
2-dimensional metric 



and use the abbreviation: — det gmn = F = + BC. The conditions for the 
negative definiteness are: 



ds"^ = dx^ + A{x)dy'^ + gmn{x)dx"^dx 



.n 



(12) 




(13) 



A>o r > 0. 



The inverse reads 




(14) 
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4 The Einstein equations 



Since we are dealing with the vacuum case, the Einstein equations reduce to 

Rij = 0, 

where Rij is the Ricci tensor of the metric (jl2|) . It is well known that, due 
to the Bianchi identities, not all the Einstein equations are independent. In 
this case it is convenient to take, as our basic equations: 

Ryy = 0, Rzz = 0, Rtt = 0, Rtz = 0. 

It is a straightforward calculation to compute the Ricci components. The 
explicit expressions read: 

Ryy = 0^ -2TA A+ri^Af -A AT= 0, (15) 

R,, = 0^2T B +2B [pf - B [2P P +C B -B C -T^ = 0, (16) 

= ^ 2r C +2C [pf - C (2P P +B C -C B -T^ = 0, (17) 

R,, = Q^2V P +2P BC - p[bc +C B -T^ = 0, (18) 

where /= ^. First of all, let us notice that eq. (fl^j) can be explicitly solved 
for A: 

where Ji and n are integration constants. Now, since A is expressed in terms 
of r, it is clear that, thanks to the identity ()20|1 . we are left with a system 
of three equations in the three unknown functions B, C and P. This system 
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looks highly nontrivial due to the nonlinearities. Nevertheless, it is possible 
to further reduce it by rewriting the equations (fTB|) . (fTTj) and (fTS|) in a more 
symmetric way. It is important to stress here that, thanks to ()16|) and ()17|) . 

is a symmetric function in the exchange of B and C. Hence, P can 
be either symmetric or antisymmetric in the exchange of B and C. In the 
following 5, C and P will be supposed to be different from zero. In fact, if 
P is zero one obtains a Kasner-like solution, while if B, or C, is zero then 
it is always possible to make a coordinate transformation in such a way that 
B ^ 0, the same holds for C. Now, it is easy to show that the system of 
equations (fTB|). (fT7|) and (fTH|) is equivalent to the following system: 

2r|H-2E-|(r-^)=0, (21) 

where Yi = B, Y2 = C, = P and = (pY + BC- It could look 



that in these equations we cannot put A= 0. However, by remembering the 
identity (f^UI) . it is obvious that one obtains A= by taking in these equations 
K = 0. Thus, we rewrote the system of the equations (fTH|) . (fT7|) and (fTH|l 
in a manifestly symmetric form: all the 1^'s obey the same equation. By 
introducing the variables rjij =Yi Yj— Yj Yi, i.e. the Wronskian of Yi and Yj, 
and supposing Vijj^ (otherwise, if, for some i and j, Vij= ^ ~ 1^- 
and the system is immediately reduced) we arrive at the following system of 
equations: 



K 



2T Vij -y^-jj = 0- (22) 
From this system, it immediately follows that 

so, after a trivial integration, it comes out that the r^jj's are all proportional: 

r]i2 ~ ??i3 ~ '723- (23) 
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As it is well known from the theory of the linear system of ordinary differential 
equations, that (|23p implies that one of the three unknown functions is a 
linear combination of the other two. It is convenient to choose P as dependent 
function. In fact, since we know that is symmetric in the exchange of B 
and C, then the only possibilities for P are: 

P = a{B±C), 

where a is an arbitrary nonzero constant. In this way, we reduced the initial 
system of four nonlinear equations in four unknown functions to a system 
of two equations in the unknown functions B and C. In general, due to the 
term ^ that couples in a non-trivial way B and C, it is not possible either 
to decouple the two equations or to further reduce the system. However, if 
one takes 

P = ±\{B-C), (24) 

then the system (j2T|l can be reduced to one equation in one unknown function. 
In fact, if eq. holds, then 

1_ ^ 1/- 



T = -{B + C)\^ = -[B + C 
Now, if we introduce 

u = B + C,v = B -C, 
then, from eq. (I2H), u and v satisfy the following two equations: 

= 0, (25) 

_ , v(---u^ = 0. (26) 

Since eq. (j^Bj) is a linear homogeneous ordinary differential equation, v can 
be expressed in a closed form in terms of u, so the system (PT|) is reduced 
to the only eq. ()25|) . Hence, once eq. (j^^ is solved for u, the other metric 
coefficients immediately follow. 
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4.1 An explicit example 

An interesting explicit example is the following: 

ds^ = dx^ + x^dy^ + ^1 + - In |a;|^ dz^ — ^1 — - In \x\j dt^ — In \x\ dz dt 



where - < ^J\x\ < e . (27) 

In this case we have: 

B + C = 2, B ~ C = \n\x\, P = -^{B -C) . 

Then it is trivial to show that u = B + C satisfies eq. (j25|) . while v = B — C 
satisfies eq. ()2fi|l for k = 4 and A follows from eq. (fTIH) by taking Ji = 0. 
It is interesting, at this point, to make a comparison with the Kasner case. 
In particular, one could ask: why the procedure to diagonalize the metric in 
the Kasner case works and in this case does not work? In this case, the first 
part of the exercise is the same as the Kasner one: Let us take the initial 
hypersurface [a; = 0] , and then without loss of generahty let 




a(0) = 1 and g„^n{0) = \ ^ ^ | . (28) 
Then we define 

dx 

and try to diagonalize fcmn(O). However, due to the non-compactness of the 
Lorentz group, the arguments of Appendix A no more apply; moreover, see 
Appendix C, one can really find examples of matrices fcmn(O) which cannot 
be brought into diagonal form. The calculation to be done is straightforward 
and will be omitted here. Then, it is clear that the differences between the 
two cases are group theoretical in nature. 
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Appendix A 



Let 



k — kaf 



^ a d e ^ 

d b f 
[e f c) 



(30) 



be any symmetric matrix, i.e. k = k""". We want to show that in the positive 
definite case, this matrix can be diagonahzed; the 3-dimensional real orthog- 
onal group is denoted by 0(3), and the superscript T denotes the transposed 
matrix. Then U G 0(3) acts continuously on k to give 



-k-U, 



(31) 



We have to show that one can always choose U G 0(3) such that the matrix 
eq. ()31|) has diagonal form. To this end we define the quantity 



J(k) = d^ + e^ + f. 



(32) 



Due to the compactness of 0(3), the minimum of J exists; we have to prove 
that this minimum leads to J = 0. Assumed, this is not the case. Without 
loss of generality we may assume that this is due to d 7^ 0, for otherwise, a 
permutation of the coordinate axes would lead to this inequality. 



Let A^ G 0(3) be defined by 

/ 

A,^ = 

V 



COS (p 

- sin I 




sm( 

COS( 





\ 


1 / 



(33) 



It holds: The inverse matrix to A,^ equals A.,^ which is nothing but A^. 
Analogously to expression (jSH), we define 



Al ■ k ■ A. 



(34) 
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and then we get with eqs. (j32j) and (|33|) up to hnear order in the Taylor 
expansion with respect to 

J{k^) = J{k) + 2{a-b) ■ d- (j). (35) 

For a 7^ 6 we are aheady finished: A small change of (p will change the value 
of J linearly with (p, so there cannot be a minimum at = 0. For a = b, this 
linear expansion does not suffice to decide, but here, the exact value is easy 
to evaluate, it reads: 

J(k^) = J(k) - 4 ■ (i^ ■ sin^ ■ cos^ (p . (36) 

There is a local maximum at = 0, and therefore, this cannot be a minimum. 
Result: the assumption d leads to a contradiction. 

Appendix B 

Equation Q, i.e., 

p + g + r=p^ + g^ + r^ = l (37) 

represents the intersection of the plane p + q + r = 1 with the unit sphere 
in the p ~ q — r— space. Thus, it must be a circle. We parametrize it by the 
angular coordinate 0. The following 3 points 

P = (1,0,0), Q = (0,1,0), i? = (0,0,1) (38) 

are obviously on this circle; in turn, this circle is uniquely determined by 
them. The center M of this circle is given by the arithmetic mean of P, Q 
and R, i.e. 

M^(iii). (39) 
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Its radius equals the distance from M to P, i.e. y2/3. So, we get the 
parametrization of eq. ^7\i as 

1 2 

p=- + -cos0 (40) 

g= ^(l-cos0) + yi73 sin0 (41) 

r = - (1 - COS0) - -y/l/S sin0 (42) 
3 

Obviously, it suffices to restrict to the 0-interval < < 27r. However, a 
permutation between the 3 numbers p, q and r can be compensated by a 
coordinate transformation (namely a related permutation of the spatial axes 
of metric 0), therefore, to get a one-to-one correspondence it proves useful 
to require additionally r < q < p. Comparing eqs. ()4H1 and (pn|) one can see 
that the inequality r < g is fulfilled for < (p < tt only. The other inequality, 
q < p further reduces this interval via the identity 

p — q = cos (j) — (sin 0) / (43) 

to 

O<0<vr/3. (44) 

Clearly, as there are six possible permutations, the length of this interval is 
2tt/6. The boundary of this interval consists of two points. The point related 
to = is the already discussed fiat space-time. The other one, related to 
= it/3, i.e. that point where p = q = 2/3, r = — 1/3 is the other axially 
symmetric solution for Bianchi type I. Here one can see what one also meets 
in other circumstances: The solutions with higher symmetry (here: axial 
symmetry) are at the boundary of the space of solutions. 
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Appendix C 



Let 



V 



( 







1 







1 



) 



(45) 



The Lorentz group 0(1, 1) is the group of all those transformations leaving 
the matrix eq. invariant. For a given symmetric matrix 



defined by 3 parameters, the one-parameter group 0(1, 1) acts continuously, 
so, from counting the degrees of freedom one could be tempted to assume, 
that one can always choose an element of 0(1, 1) such that D becomes zero. 
However, this is not the case. For our purposes it suffices to give an example: 
For 



put into eq. pUj) . no diagonalization is possible. 

Let us look from another side: The trace of k, namely the expression^ 
A—B, and the determinant, namely AB — D^, are invariants of it with respect 
to 0(1, l)-actions. Again, the counting is misleading: 3 free parameters in 
fl46|) . a one-parameter gauge group, so these two invariants should suffice for 
an invariant characterization. But this is not true: both A = B = D = 1/2 
and A = B = D = lead to a vanishing of both invariants, whereas no 
element of 0(1, 1) can be given that transforms the one into the other. This is 
analogous to the discussion in jH]: two objects are different, but no continuous 
invariant exists to distinguish between them. Here it holds: every continuous 
invariant of k can be written as a function of trace and determinant only. 
^The minus sign in front of B is due to the minus 1 in eq. H45|) . 




(46) 



A = B = D = 1/2 



(47) 
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Appendix D 



Let a metric be given as 

ds"^ = dx^ + dy"^ + 2dzdt + a{x) df (48) 

where a{x) is any free function. Due to the off-diagonal term dz dt, a{x) may 
have zeroes without leading to a singularity there. 

We denote (x, y, z, t) hj x\ i = 1, . . . 4. The only component of the Ricci 
tensor Rij which does not vanish identically, is 

1 / d'^a 1 da\ , . 

The only components of the Riemann tensor Rijki which do not vanish 
identically, are 

1 d'^a X da , . 

^'''' = '2'd^^' ^'''' = ~2'd^- 

This statement is meant, of course, only "up the usual antisymmetries". 

As a result of eq. (j3n|) we find: Metric (jlH|) is flat if and only if the 
function a is a constant. 

To find out all non-fiat vacuum solutions of the Einstein field equation of 
the form (pSj) . one has therefore to solve -R44 = using eq. with a non- 
constant a{x). The result is, after a possible redefinition of the coordinates 
t and z, be expressible as 

a{x) = c ± Ina; (51) 
where c is a given constant of integration. 



Let us calculate the curvature invariants of metric (j48|) : Let I be any 
polynomial invariant like K^^ Rij . Then / depends on the one coordinate x 
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only. To calculate one special value I{x) we make the following construction: 
We replace, for any positive real e, the coordinate thyet and the coordinate 
z by z/e. This does not change the form of metric (|lHj) . only the function 
a(x) is now replaced by • a{x). 

In the limit e — 0, we meet the flat spacetime having / = 0. But I{x) is 
a continuous function, and as invariant it does not change with e, therefore: 
Every polynomial curvature invariant for metric PHj) identically vanishes. 

Note added 

The paper [9] and ours are different in scope, but the discussed metrics have 
much overlap: [9] presents the most general metric that depends on just one 
coordinate and cannot be diagonalized. The metric is a generalization of the 
Levi-Civita, or Kasner metrics. The authors of [9] also analyzed the global 
structure of the spacetime described by this metric - it has closed timelike 
curves, without a Cauchy horizon, and the question of whether such metric 
can represent the "exterior" of some "tube of matter" , the answer is that in 
general an energy condition is violated. 

In [10], several results on diagonalization procedures can be found, too. 
According to M. MacCallum, the solution in the abstract of the present paper 
is just a special pp-wave and appears as eq. (22.5) in the new edition of ref. 
[5] (with a = 1 and p = Inx). It also represents a special case of solutions 
already given in refs. [1] and [4]. In refs. [11] and [12], similar solutions have 
been discussed, too. 
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